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Abstract 

The Blandford-Znajek (BZ) mechanism describes a physical process for the energy extraction 
from a spinning black hole (BH), which is believed to power a great variety of astrophysical sources, 
such as active galactic nuclei (AGNs) and Gamma ray bursts (GRBs). The only known analytic 
solution to the BZ mechanism is a split monopole perturbation solution up to O(a^), where a is the 
spin parameter of a Kerr black hole. In this paper, we extend the monopole solution to higher order 
~ 0(A). We carefully investigate the structure of the BH magnetosphere, including the angular 
velocity of magnetic field lines D, the toroidal magnetic field B^’ as well as the poloidal electric 
current I. In addition, the relevant energy extraction rate E and the stability of this high-order 
monopole perturbation solution are also examined. 
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I. INTRODUCTION 


Within the framework of force-free electrodynamics, Blandford & Znajek (1977) investi¬ 
gated a steady-state axisymmetric magnetsphere surrounding a spinning black hole and put 
forward that the rotation energy of a Kerr black hole could be extracted in the form of Poynt- 
ing flux via magnetic helds penetrating the central black hole ([H |2])- General relativistic 
magnetodynamics (GRMD) simulations of split monopole magnetic helds ([ 3111 ]) show that 
the analytic monopole perturbation solution makes good matches with the numerical simula¬ 
tions, especially for slowly rotating black holes. General relativistic magnetohydrodynamics 
(GRMHD) simulations ([5l|6|) indicate that, in the polar region, the monopole perturbation 
solution gives a good description of the magnetic held conhguration as well as the angular 
distribution of energy how, even when black holes rotate mildly rapidly. However, the mono- 
pie solution [T] is accurate only up to O(a^), where a is the black hole spin parameter. For 
even more rapidly rotating black holes, higher order perturbation solutions are of greater 
astrophysical interests. Tanabe and Nagataki [7] extended the monopole perturbation so¬ 
lution to the order of 0(0“^). Their solution gave a better approximation to the numerical 
simulation. Unfortunately, they mentioned that their results are not fully self-consistent, 
since their perturbation method breaks down at large distance from the central black hole. 
Hence, it is necessary to hnd self-consistent higher order perturbation solutions to the BZ 
mechanism. 

To get self-consistent solutions, we need to solve a nonlinear second-order partial diher- 
ential equation, which requires two boundary conditions. It should be noted that boundary 
conditions to be imposed are still not well understood ([THTT]). Blandford and Znajek jl] 
imposed the Znajek regularity condition (ca) on the horizon as the hrst boundary condi¬ 
tion. The second one requires that the perturbation solution should match the asymptotic 
solution in the flat spacetime at inhnity ([I3|). Unfortunately, the second boundary condi¬ 
tion is usually unavailable when investigating higher order perturbation solutions. Recently, 
Pan and Yu [H] proposed that the physical constraint, i.e., solutions should be convergent 
from the horizon to inhnity, could be exploited as the second boundary condition. With 
the Znajek horizon regularity condition and this new convergence constraint, perturbation 
solutions could be uniquely determined. Following the approach of Pan and Yu HI , we ex¬ 
tend the monopole perturbation solution to the order of O(a^). Note that the perturbation 
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method we adopt is different from [7]. Our method works well at any distance from the 
central black hole. 

Some earlier analytic works (dSHHl) concerned the stability of jets launched by the BZ 
mechanism because of the screw instability of the magnetic field. However, such instability 
was not found in recent simulations (e.g. [IHH2Q]). The possible reason for the discrepancy 
is that the Krustkal-Shafranov (KS) criteria is used in these works, without taking account 
of the stabilizing effect induced by the magnetic field rotation (|THl l2ll 122] )• With the 
high-order perturbation solution obtained in this paper, we also briefly study the stability 
of the split monopole perturbation solution of the order of 0{a^), taking the magnetic field 
rotation into consideration. 

The paper is organized as follows: basic equations governing stationary axisymmetric 
force-free fields around Kerr black holes are described in section II. We discuss the pertur¬ 
bation solutions of second-order and fourth-order obtained by our newly proposed method 
in section III. Summary and discussion are given in section IV. 


II. STATIONARY AXISYMMETRIC FORCE-FREE FIELDS AROUND KERR 
BLACK HOLES 


In this section, we briefly recap basic equations governing stationary axisymmetric force- 

free fields around Kerr black holes (see [13] and references therein for more details). We 

adopt the Kerr-Schild coordinate (e.g., McKinney and Gammie [5]) with the line element 

, 2 2r\ ,2 /4r\ , , 2r\ ,2 ^,^2 4arsin^6' ,, , 

ds^ = — 1 1 — — I -I- 1 — I drdt -|- (1 -|- — j dr^ + YldO^ -^- dcjydt 


E/ 


E/ 


2a ^ sin^ Odcpdr -|- sin^ 9 


A + 


2r(r^ -|- 


d4P‘ , 


( 1 ) 


where E = cos^ d, A = — 2r -|- a^, and = E sin d . 

The energy momentum tensor for the force-free field is dominated by the electro-magnetic 
field, which can be written as = T,((atter + ~ ^em = ^ — \5^''F°'^Fa^^ where 

the antisymmetric Faraday tensor is defined as F^^ = and A is the 4—potential 

of electromagnetic field. We define the angular velocity of the magnetic field VL[r,d) as 
follows. 
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dAf A 
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It is evident that = 0 for the axisymmetric and steady state force-free held. The non-zero 
parts of the Faraday tensor are listed below: 


r(p 


tr 


rt 



( 3 ) 

1 ^t9 ^9t ^-^(p,9 1 

( 4 ) 

= -For = . 

( 5 ) 


Note that the force-free held is specihed by three quantities, i.e., A^{r,9), and 

B'^{r,9). Once they are specihed, the force-free held is uniquely determined. 

Note that T® = and T\ = The energy and angular momentum con¬ 
servation equations = 0 and = 0 can be cast as and 

= {y/—gF^^)^ 0 A^^r, respectively. These two equations indicate that O 
and are functions of A^, viz. 


O = n{A^) , v^F^’' = liA^) , 


( 6 ) 


where the angular velocity of magnetic held O and the poloidal electric current / are to be 
specihed. Substituting Equations ([^, @, (|^ and § into the equation F®'" = 


we can readily arrive at 


= 


IS -|- (20r — a) sin9A^fi 


( 7 ) 


AE sin^ 9 

This is an important relation that connects the toroidal magnetic held B^ with the functions 


A^{r,9), and 

The remaining momentum conservation equations in the r and 9 direction = 0 and 
= 0 are actually equivalent and read 




+ FreI\A^) + J = 0 , (8) 


where the prime denotes derivative with respect to A^. The three functions A^{r, 9), 11(^0), 
and I{A^) are related by the above nonlinear equation (|^, which is also widely known as 
the Grad-Shafranov (GS) equation ([HIES]). 


III. FOURTH-ORDER PERTURBATION SOLUTIONS 

Since the Farady tensor depends on the first order derivative of it is clear 
that the GS equation @ is actually a second-order partial differential equation 
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for Aip. The solution can be attained when complemented with two boundary conditions, 
i.e., the Znajek horizon regularity condition ([12]) and the convergence constraint (PEI)- 
The zeroth-order monopole solution can be readily obtained when the black hole is non¬ 
rotating, i.e., a = 0. When the spin parameter a 7 ^ 0, we expand the GS equation in 
terms of a. To get the second-order perturbation solutions, we ignore all terms in the GS 
equation that are higher than the order of O(a^). Based on the second-order solutions, the 
fourth-order perturbation solution can be achieved in a similar way. 

The zeroth-order monopole solution around non-rotating black hole can be explicitly 
written as (P), 

Go = 0, Bq =0, = Aq = — cos 9. (9) 

Since G and y/—gF^^ are functions of A^, we can expand them, accurate to the order of 
0 (a'^), as 

G = G(74^) = auJi{A^) = qui^Aq A 2 ) -\- a^u^^Ao a‘^A2), 

= I{A(j)) = aii{A(j)) a^is[Aff,) = aii(Ao -|- o?'A 2 ) -\- a^is{Ao a^y42), (10) 

where G, cui, 0 J 3 , /, G, G are unknown functions of A^j, to be specified self-consistently. The 
entire fourth-order perturbation solutions can be expressed in a more compact form as, 

A^ = Aq CL^ A 2 A 4 O(a^) , 

G = aGi -|- a^Gs -|- 0{a^) , 

^F^'^ = ah + a^h + 0{a^) , 

= aBi + a^Bs + 0{a^) . (11) 

It should be noted that G„ and Un, h and G {n = 1,3) are related by 

Gi = Ci;i(Ao) , G 3 = uj'~^[Aq)A2 + Ci;3(Ao) ) 
h = h{^o) 1 h = *i(^o)^2 + G(^o) ) (12) 

where the prime designates the derivative with respect to Ao- 

A. Second-order Perturbation Solutions 

We can get the second-order perturbation solutions by expanding the GS equation ([^ 
to the order of O(a^). It is interesting that the original BZ monopole perturbation solution 
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could be naturally achieved with our convergence constraint. Expanding Eq.([^ to the order 
of 0 {a^), we have that 

r^/i = sin6*y4o,0(l — 2rf2i) — sin^ 9Bi{r‘^ — 2r)r^. (13) 


According to the Znajek horizon condition ([12]), the toroidal held Bi should be well-behaved 
on the horizon (r = 2), then r = 2 must be a root to equation = sin^ 6^(1 — 2rfli). So 
we have 


ii = Ii = sin^ 6^ ( - — 111) , 


--(5 (i - ^ 


The GS equation ([^, accurate to the order of O(a^), can then be cast as 

CA2 = S{r,9) , 


(14) 


(15) 


where the operator 


and the source 


1 a / 2\ a 1 d 1 d 

sin 0 c?r V r ) dr 89 sin 9 89 


S{r,9) = 4sin0cos0 ( ^ “ 2sin0cos0^ ^ 

1 1 




(16) 


(17) 


According to Blandford and Znajek [T], the condition for the existence of solution is that 
the following integral, 

J^(a0)I 


/ dr d9- 
h Jo r 


(18) 


should be convergent. The convergence condition requires all the terms in S{r,9) of the 
order of 0(1) should vanish, i.e., 

1 ' 


0 = 4 sin 6* cos 6* fli — - M- sin g . 


(19) 


Consequently, we have 


fii — (jj\ — — , 


Ii=ii = - sin^ 9 , 

O 


V8 2r 


( 20 ) 
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It is interesting to note that all physical quantities of the order 0(a) are already obtained 
before we actually solve the complicated GS equation. The second-order part of A^, i.e., 
A 2 , can be obtained by the following equation, 

£^2 =- 2 sin 6 'cos 6 '^ . ( 21 ) 

V2r 

It is straightforward while tedious to check that this equation has the following variable 
separable solution(|l]) 

A 2 = R(r) sin^ 9 cos 9, (22) 


where 




12 


r\ 11 1 r 

^ 2~ Y 


and 


+ 


LUix) = — dt 


Lio - - In 1 - - In hr 


r^(2r — 3) 


ln(l — tx) 


t 


(23) 


(24) 


The value of the function R(r) at the horizon, r = 2, is of particular importance. Explicitly, 
it is 


Rr =2 = R(r = 2 ) = 


bvr^ — 49 
72 


(25) 


B. Fourth-order Perturbation Solutions 


Once the second-order perturbation solutions are known, the fourth-order perturbation 
solutions could be obtained by further expanding the GS equation to the order of 0(a"^) . 
Accurate to 0{a'^), Eq.([^ is 

r ^/3 -|- cos^ 9Ii = sin 6 * [A 2 ,e(l — 2 rGi) — AQ^Q2rVL^ 

— sih 6* B^(r‘^ — 2r)r‘^ + Bi(r‘^ — 2r) cos^ 9 + Bir'^ . (26) 


The toroidal held B^ should be well behaved on the horizon. Subsequently, we can get 

h -h sih 9u:^ = ^ 

and the toroidal held B^ is 


(27) 


sin^ OB, = 


— 2r 


2 sin^ 9 


-COs + 


sih0/l 1 
8 ^ ^ 


1 \ . sin 0 cos^ 9 ' 

Sint* A 2 ft - o 


4r 


hos 9 


A2 + h 


(28) 
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where we have made use of Eq.(12) and (20). The GS equation (|^ of the order of O(a^) is 
of the following form, 


— UJ\ 




^3 


,r 


1 + - sin^ 6 * + 2 r sin^ 9Bi 

rj 8 


-(Ui 


• sin 0 ( 1 + ^ ^ ( ^^ 2,0 + 1^3 sin 9^ + 


rJ \8 


sin^ 9 cos^ 9 


J ,0 


+ 2 r sin^ 93^ 


(29) 


J , 6 » 


+r^ sin 9Bi 
( sin 6 ' 


-^2 + * 3 (^ 0 ) + T cos 6 * sin 0 (^cos^ 9Bi + r^B^ 


2 cos^ 9 \ / . 2 /I n cos^ 9 cos"^ 9 \ 

— 1 —I ^2,r H-3 . n ^2,r + SlU 9B'i + . „^2,0-^ 

y 4r r-^smt7 J ^ y r^sint^ r° J ^ 

The convergence condition requires all source terms of the order 0 ( 1 ) should vanish, i.e., 

0 = a;3 sin^ 6*^ + a;i(sin^ 6*ci;3)_6» +sin 6*51*3(740) + 7 sin 0 cos 6*r^i?3 . ( 30 ) 

Vo /,0 ’ 4 

The above equation could be further simplihed as 

COS 0 1 

a;3(sin^6'),6i + (sin^ ^ws),^ = *3,0 + = ^^(sin^ 6'*3),0 <t 4 sin^ 6*a;3 = *3, ( 31 ) 

sm 9 sm 9 


where we have used the result of Eq.(20) and (28). Together with Eq.(27), we have that 

1 / Rr=2 1 


*3 = 


8 


H-) sin'^^ H-sin^ 6 *, 

32/ 32 


1 / Rr=2 1 \ • 2 /I ^ ^ 

^ 2 V 8 32/ 32 32 


(32) 


With the help of Eq.(12) and (20), we finally arrive at 

a 


n — V2(74<^) — o + Q (^3) 




(33) 


IV. DISCUSSION AND SUMMARY 


A. Discussion 

The angular distribution of the fourth-order angular velocity O and poloidal eletric current 
I on the horizon is shown in Fig|^ For comparison, the corresponding simulation results 
are also available (cf., Fig.l and 2 of [3]). Both simulations and our analytic solution imply 
that O = Oh /2 is a rather good approximation for a wide range of black hole spins (at 
least for a < 0.9), where Oh = a/2(l -|- \/l — a^) is the angular velocity of the central BH. 




























FIG. 1. Angular distribution of the ratio Q/^Ih and the electric current I on horizon (r = 2), 
where we keep the angular velocity of the BH accurate to fourth order, i.e., Qh = a/4 + a^/16 . 

The fourth-order poloidal electric current I also shows better agreement with the simulation 
result than the second-order one, especially for large spins. 

The energy extraction rate, which is defined as E = —27i Jq y/^T\d9 = 27 t f I{A^)Q{A^)dA^ 
(Pill El), could be written as 

E = 27ra^ J iiUidAo + 27ra^ J + uAs + {iiUJiA2ydAo 

= 2710? J iiUidAo + 2770? J g (^3 + 9u:?)dAQ 

71 9 (56 — 371^)71 4 , . 

= —o^ + ^- —o^ 34 

24 1080 ^ ^ 

where the prime denotes derivative with respect to Aq. Note that the second term on the 

right hand side only depends on the combination, Z 3 -f-sin^ which can be specified by the 


Znajek horizon condition (i.e., Eq.(27)). In fact, this coincidence explains why Tanabe and 


Nagataki [7] could obtain the correct energy extraction rate E without explicitly solving G 
and J. 

The stability is another interesting issue. Some analytic works ( [TSHTT] ) implied that 
the screw instability may occur in the monopole perturbation solution due to the Kruskal- 
Shafranov criterion. But no instability was noticed in time-dependent GRMD (e.g. Pi) 
or GRMHD simulations (e.g. [HEn]). To understand the discrepancy between analytic 
and numerical works, Narayan et al. [25] (and [HEI]) pointed out that Kruskal-Shafranov 
criterion may not be appropriate for jet stability analysis, since it neglects the stabilizing 
effect of the rotation of magnetic held lines. According to the analysis of Tomimatsu et al. 
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[21], which takes the held rotation into account, the monopole perturbation solution is 
possibly unstable only when Our fourth order solution (i.e., Eq.(33)) means 


that 


O > 





(35) 


Obviously, the fourth-order monopole perturbation solution is stable and is consistent with 
numerical simulations. 


B. Summary 


Two major difficulties are encountered in solving the GS equation (|^: 1) it is a highly 
nonlinear second order partial differential equation; 2) two proper boundary conditions are 
necessary to uniquely specify the solution. The nonlinearity could be partially removed 
by the perturbation technique. To hx the boundary conditions problem, we impose the 
regularity condition on the horizon (Eq.Q) and the convergence constraint (Eq.(18)). The 
latter one actually serves as the boundary condition at inhnity. With these two boundary 
conditions, we re-establish the split monopole solution to the order of O(a^) and get the new 
perturbation solution up to the order of 0(a‘^). By taking account of the stabilizing effect 
of held rotation, we prove that the fourth-order monopole perturbation solution is stable 
against the screw instability. 
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